Distance-redshift from an optical metric that includes absorption 
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We show that it is possible to equate the intensity reduction of a light wave caused by weak absorp- 
tion with a geometrical reduction in intensity caused by a "transverse" conformal transformation 
of the spacetime metric in which the wave travels. We are consequently able to modify Gordon's 
optical metric to account for electromagnetic properties of ponderable material whose properties 
include both refraction and absorption. Unlike refraction alone however, including absorption re- 
quires a modification of the optical metric that depends on the eikonal of the wave itself. We derive 
the distance-redshift relation from the modified optical metric for Friedman-Lemaitre-Robertson- 
Walker spacetimes whose cosmic fluid has associated refraction and absorption coefficients. We then 
fit the current supernovae data and provide an alternate explanation (other than dark energy) of 
the apparent acceleration of the universe. 
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I. INTRODUCTION 

The concept of an optical metric was introduced by Gordon [l| who proved that solutions to Maxwell's equations 
in a curved spacetime filled with a fluid whose electromagnetic properties are described by a real permittivity e(x) 
and a real permeability /i(cc) (and consequently by a real refraction index n(x) = y^e/i) could be found by solving a 
modified version of Maxwell's equations in a related optical spacetime with vacuum values for the permittivity and 
permeability, i.e., with e(x) = 1 and /j,(x) = 1 [see Eq. ([64]) ]. There is a direct relation between solutions in these 
two spacetimes, the physical spacetime with an index of refraction and its physical metric, and the optical spacetime 
with n — 1 and its optical-metric. Gordon's original metric accounted for refraction only. An open question was 
whether or not absorption could also be accounted for by modifying the spacetime geometry. In [2| we made such 
a modification and showed that Gordon's optical metric could be generalized to include absorption by allowing the 
metric to become complex. In Q we defined a complex refraction index N = n + in, and distinguished two different 
cases, i.e., strong and week absorption. In the case of strong absorption the real and imaginary parts of the index of 
refraction are of the same order (k ~ n), the eikonal S(x a ) has an imaginary part, and the optical metric is complex. 
In the weak absorption case (I; « n), the eikonal S(x a ) can be taken as real, and the real part of the optical metric 
(Gordon's original metric) remains as the significant geometrical structure. In this paper we show that for weak 
absorption another generalization of Gordon's metric exists that can account for absorption as well as refraction. 
The generalization amounts to a transverse conformal transformation (see Sec.[llJ) of Gordon's original optical metric. 
This new optical metric remains real but does include absorption. Including absorption via a transverse conformal 
transformation requires an optical metric that depends on the eikonal of the wave itself. 

The contents of this paper are as follows: In the next section we define transverse conformal transformations. In 
Sec. lIIII we consider the geometrical optics approximation and relate wave properties in the two spacetimes, the physical 
and optical. In Sec. lIVI we compare the effects of the transverse conformal transformation with absorption and relate 
the conformal factor a to optical depth r. In Sec.[V] we construct the optical metric that accounts for absorption in 
Friedman-Lemaitre- Robertson- Walker (FLRW) spacetimes. In Sec. lVIl we include refraction by generalizing Gordon's 
metric [l[ to including absorption. In Sec. lVIII we derive the refraction/absorption corrected distance-redshift relation 
from the generalized Gordon's optical metric in FLRW spacetimes. In Sec. lVIIll we fit the current supernovae data 
with the Hubble curve of a cosmological model in which the cosmic fluid has both an associated refraction index n and 
a constant opacity parameter a, and provide an alternate interpretation of the apparent acceleration of the cosmic 
expansion. In Sec. lIXI we give our conclusions. 
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II. THE "TRANSVERSE" CONFORMAL TRANSFORMATION 

A conformal transformation Q is defined as a rescaling of the metric and is usually written in a form similar to 

9ab = e 2a g ab , (1) 
where a(x a ) is an arbitrary scalar function defined on the spacetime manifold. The vacuum Maxwell equations 1 

F[ab,c] = 0, 

V b F ab = 0, (2) 
are formally invariant under a conformal transformation, i.e., 

F[ab,c] = 0, 

i/ b F ab = 0, (3) 
provided that the covariant electromagnetic field tensor F a b transforms as 

F ab = F ab , (4) 



and hence the contravariant field transforms as 



pab = e -^ p ab_ ( 5 ) 



The above is a purely mathematical observation and simply says that if you have a solution to Maxwell's vacuum 
equations in one spacetime then you have a related solution to Maxwell's vacuum equations in any conformally related 
spacetime. Nothing is being said by this about the existence of any new physically relevant field. One can and does 
make good use of such transformations to compactify spacetimes and to analyze fields at "oo", see e.g., [4|. However, 
by modifying these conformal transformations to be "transverse" [see Eqs. (TlQ)) and (flTT) ] we are able to construct new 
and useful solutions to Maxwell's equations. 

Because the two metrics are defined on the same differentiable manifold we can make unique correspondences 
between events, world lines, and various fields in these two spacetimes. We call the original metric and manifold 
physical spacetime and the manifold with the conformally transformed metric the conformal spacetime. For example 
a normalized fluid 4- velocity, u a (x), defined in the conformal spacetime would be related to the corresponding physical 
fluid 4-velocity by u a (x) = e^ cr u a (x). If we compare observed properties of radiation fields that are related as in 
Eqs. (|4]) and ([5]) and as seen in these two spacetimes by corresponding observers, we find that energy fluxes are 
diminished by a factor of e~ 4<T . This net effect can be looked at as the result of the area of the radiation's wave front 
changing by a factor e 2<T , and the energies and rates each changing by a factor e~ CT . We will make use of this intensity 
reduction to account for absorption in physical spacetime. However, because a pure conformal transformation has the 
undesirable property of altering a wave's frequency and wavelength it must be appropriately modified to represent 
absorption, which obviously does not. 

We assume that our physical spacetime is filled with two vector fields. The first is the normalized velocity field 
u a (u a u a = —I) of a fluid whose linear and isotropic optical properties we know and the second is a null field k a 
(k a k a = 0) corresponding to the eikonal of a given electromagnetic wave (see Sec. lIIII for details of the geometrical 
optics approximation). These two fields define a 2-dimensional timelike subspace at each point of the manifold. We are 
hence able to decompose the tangent space at each point of the 4-dimensional spacetime manifold into two orthogonal 
2-dimensional subspaces, one timelike g\\ a b and spanned by the pair (u a ,k a ) and the other spacelike g± a b- The full 
metric can be decomposed in terms of these two orthogonal pieces as 

9ab = 9\\ab + gjLab- (6) 

We can give a simple expression for g^ ab by defining a second null vector l a lying in the two dimensional timelike 
subspace spanned by it° and k a as 

h a ii a 

(7) 



2(u-k) 2 (u-k)' 



1 Square bracket [■] or parenthesis (■) symbolize respectively complete anti-symmetrization or symmetrization of the enclosed indices, and 
V with a subscript is covariant differentiation. 
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from which it follows that 



ti a = 0, tu a = — — — — TT", tk a = 1. (8) 

2{u ■ K) 



The metrics on the orthogonal 2-dimensional surfaces can be written as 



9\\ab = 2£( fc b ), 

9±ab = 9ab-2i(ah)- (9) 



We conformally transform the 2-dimensional spacelike subspace only and arrive at the desired optical metric 

9ab = g\\ab + e 2lJ g±ab, 

T gab + (l-e 2 °)2e (a k bh (10) 



with inverse 



g ab = gf + e-^gf, 



= e- 2a g ab + {l-e- 2a )2£( a k b \ (11) 
By straightforward tensor algebra we find that 

det g = e 4CT det g, (12) 

We would have obtained det g = e 8tT det g were we conformally transforming the entire metric g a b as in Eq. ([T|) . We call 
these transformations transverse-conformal because they scale directions transverse to a wave's propagation direction 
k a as seen by an observer u a moving with the optical fluid. We find that the vacuum Maxwell equations ([2|) remain 
invariant in form under the transverse conformal transformation of Eq. (|10[) . i.e., 

F[ab,c] = 0, 

V b F ab = 0, (13) 
provided that the covariant electromagnetic field tensor F a b transforms as 

F ab = F ab , (14) 
and that the contravariant field tensor F ab defined by 

pab _ ~ g ac~bdp cd (15) 

is constrained to satisfy 

pab = e -2a p ab ^g) 

This constraint requires that two of the following three terms, i.e., F§ b and F£ b , originating from Eq. (fT5|) vanish: 

pab = pab + p -2apab + pab ^ ^ 



where 



F o6 = 2t [a (F cd k c t d ) , 
Ff = -2k [a F b] J c - 2£ [a F b] c k c - 2F afc , 

pab ^ pab _ pab _ pab _ ^ 



Consequently F a b must satisfy 



and 



F cd k c £ d = 0, (19) 
F ab = -2k [a F b] J c - 2£ [a F b] c k c . (20) 
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As we will see in the next section [see Eq. (|32[) ] a radiation field whose eikonal generates the null field k a satisfies 
these constraints. For a familiar example, choose the physical metric to be Minkowskian, i.e., 

ds 2 = -c 2 dt 2 + dr 2 + r 2 (d6 2 + sin 2 6 d(/) 2 ), (21) 

with a stationary optical fluid u a = (1/c, 0, 0, 0), and radial null geodesies k a = k (f/c, 1, 0, 0). From Eq. ([7]) we find 

r = i (4,1,0,0). (22) 

Equations (|T9)l and (f2TJ)) then give 

Fqi = —Fio — 0, ^23 = —F32 = 0, (23) 

i.e., F ab is a transverse field, hence motivating the designation of Eq. (TIT)]) as a "transverse" conformal transformation. 
Equations (|19[) and (|20[) simply express the transversality condition for an arbitrary wave. In the following section 
we review the geometrical optics approximation and relate electromagnetic quantities in the physical and optical 
spacetimes. 



III. GEOMETRICAL OPTICS APPROXIMATION 

Following @, S Q we write the covariant (and metric independent) field tensor as 

F ab = Re \e lS ' x ° ( A ab + ^B ab + O(X 2 )) \ , (24) 



1 

where Ao is a wavelength related parameter, S(x a ) is the so-called eikonal function and is real, and Re{-} stands 
for the real part. The A ab term represents the geometrical optics approximation and the B ab term is its first order 
correction. Defining the unitless (also metric independent) wave vector k a — d a S and inserting Eq. (|24[) into the 
vacuum Maxwell equations we obtain to order Xq 1 

A[ ab k c ] = 0, 

A ab k b = 0, (25) 

and to order Ao 

d[ a A bc ] + fc[ Q B hc ] = 0, 
V b A ab + B ab k b = 0. (26) 

Equations (|25[) tell us that k a = g ab k b is tangent to null geodesies of g ab , i.e., 

k k a = 0, 

k b V b k a = 0, (27) 

and that A ab is of the form: 

A ab = -2k [a £ b] , (28) 

where £ a is spacelike and constrained by E a k a = with the remaining gauge freedom (to order Ao°) £ a — * £a + f(x)k a . 
Since k a £ a = 1, we can use this freedom to choose £ a such that £ a £ a — or equivalently £ a u a — 0. For this choice, 
£ a is (up to a factor lu^ 1 ) the amplitude of the electric field seen by observers at rest with respect to the fluid u a . 
It is the geometrical optics approximation, i.e.,Eq. (j28|) that makes the Maxwell field of Eq. ([24]) satisfy the needed 
transverse constraints of Eqs. (| 1 9|1 and (|20|1 and hence allows us to introduce the transverse conformal transformation 
of Eq. (fT0| . Before doing so we first finish the geometrical optics approximation for F ab in the physical spacetime. 
Equation (f2"6")l tells us that the order Ao 1 correction to geometrical optics is of the form 

B ab = 2(£ M] - k [a V b] ), (29) 

with a remaining gauge freedom, to 0(Xq 1 ), T> a — * T> a + g(x)k a and also gives the propagation equation for £ a 

£ a + e£ a = ^-(V b £ b + k b V b ), (30) 
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where '-'is the affine parameter rate of change along the null ray, £ a = fc b Vb£ a , and 9 is the expansion rate of the 
null congruence, 9 = V a fc a /2. By splitting £ a — £e a into a scalar amplitude £ and a unit polarization vector e a e a = 1 
(where * is complex conjugation), the transport equation for the amplitude £ becomes 

£ + £9 = 0. (31) 

The geometrical optics approximation, i.e., the 0(Ao°) term in Eq. (j24|) , hence simplifies to 

F ab = -2Re{£e zS/x °k^ a e b] }. (32) 

As we have indicated above, because both the physical metric g a b and optical metric g a b of Eq. (|10[) are defined on the 
same manifold we can compare properties of a common field such as F a b in both spacetimes. For the above geometrical 
optics field all covariant quantities such as F a b,k a ,£ a , and £ a in the optical spacetime are exactly the same as 
F a b, k ai t a , and £ a in the physical spacetime. All contravariant components that lie in the k-£ plane are also unchanged, 
i.e.,fc° = k a ,£ a = £ a , and u a = u a . Consequently, quantities such as affine parameters, frequencies and wavelengths 
of the waves are the same in both spacetimes. However, transverse contravariant components, i.e., components in the 
orthogonal 2-dimensional spacelike surface are scaled by the conformal factor e~ 2rT , e.g.,£ a = e~ 2a £ a . This changes 
the magnitude to £ = e~ a £ with a new unit polarization vector e a = e~ a e a and the expansion parameter 9 of Eq. (|3"Tj) 

to e = e + &. 

The time averaged 4-flux seen by an observer moving with the optical fluid is in general 

S a = ^-Re ^H ac F cb - - A 5%H dc F cd ^ u b . (33) 

In the physical and optical spaces they are related by 

S a = e- 2a S a = -e- 2a — {u b k b )££k a , (34) 

and from the 3-D Poynting vector S 1 " = (g ab + u a u b )Sb we find magnitudes related by 

S± = e- 2rJ S ± = e- 2a ^-££{u b k b f. (35) 

This result says that the intensity of a monochromatic wave can be reduced by a factor e~ 2a at any point in spacetime 
by a transverse conformal transformation without altering the wave's frequency or wavelength. In the next section 
we equate this reduction of energy flux with absorption. In Sec. I VII when we combine absorption with refraction we 
reverse the process and assume that physical spacetime possesses the absorbing material and the transverse conformal 
transformation to the optical spacetime removes it. 



IV. A TRANSVERSE CONFORMAL TRANSFORMATION VERSUS ABSORPTION 

The attenuation coefficient k u (cm 2 • g" 1 ) is defined by looking at the amount of energy absorbed, dE Vl in time dt 
by a pencil beam of radiation as it passes through a small slab of material of density p (g • cm -3 ), cross-section dA 
and length dl (see e.g. Q) 

dE u = {pK v )I v dQ,dvdtdAdl. (36) 

All quantities are defined in the local co-moving frame of the optical fluid u a . The specific intensity 1^ (xdSldv) 
measures the wave's intensity directed into solid angle dQ and within the frequency range v to v + dv. In the absence 
of emission the specific intensity I v (t), after traveling an optical depth r = J pK v dl along a "characteristic" direction, 
differs from the value J„(0) it would have without absorption by 

I V ( T ) = J„(0)e- T . (37) 

Evaluating r can be complicated because the frequency in the integrand is continually Doppler shifted due to the 
non-static nature of the optical fluid and/or the curvature of spacetime. The specific flux S u at any point is the first 
moment of specific intensity l Vl i.e., 

S v = 2ti J fj,I v dfjL. (38) 
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For any I v that has a delta function dependence on direction, e.g., a geometrical optics wave, the corresponding value 
of the flux S v (r) seen by the observer in the presence of absorption is similarly related to the non-absorption value, 
i.e., 

S v (t) = S u {0)e- T . (39) 

Comparing Eqs. (|35|) and |39|) we find that a conformal factor a reduces the flux by the exact same amount as 
absorption if 

a = \t. (40) 

If this conformal factor is to be unique, then either a single frequency is present at any spacetime point or the 
attenuation is "grey", i.e., the opacity a = pK v is frequency independent. Even if there is only a single frequency 
present the frequency dependence in the absorption makes the calculation complicated. We must follow the frequency 
of each wave, starting from the source, as it is red/blue shifted and selectively absorbed while moving through the 
optical medium. In general a different conformal factor would exist for each source frequency. However, for the "grey" 
case the a(x a ) is unique. For a geometrical optics wave emanating from the world line of a point source, a(x a ) is 
defined on forward light cones (surfaces on which the eikonal S remains constant), i.e., 

a(x a ) = ±J X p[x ,b ) Kv (x ,h )dl. (41) 

The integration is performed along the null geodesic connecting the emitting event x® and the spacetime point x a . 
The density, attenuation coefficient, and spatial length element, respectively p, K v , and dl, are measured in a sequence 
of local inertial frames (u a oc 5§) which are at rest with respect to the material along the null geodesic. We caution 
the reader that a(x a ) might not be globally defined, or even uniquely defined. The integral above is only defined for 
x a which lie on forward light cones of the source. If the light source is turned on/off at some time (star birth/death) 
a will only be defined on part of the spacetime manifold. However, its value can be taken as zero on the remainder. 
Furthermore, we might have multiple null geodesies connecting the emitting event S and receiving event O. For an 
example, when an Einstein ring occurs in a gravitational lensing configuration we have, in principle, an infinite number 
of null geodesies connecting S and O. These rays might pass through regions with different p's and k^'s and therefore 
give us different er's. Rather than limiting the domain over which the conformal transformation is defined to keep 
it single valued, we can extend the manifold to multiple layers and make the conformal factor unique on each layer. 
This is in direct analogy with Riemann's extension of the domain of a multi-valued complex function on R 2 to a 
Riemann surface on which its value is unique. This extended domain could be truly convoluted as for example in an 
inhomogeneous cosmology where strong lensing is prevalent. 

The simplest example to illustrate a transverse conformal transformation is a plane electromagnetic wave traveling 
in Minkowski spacetime filled with stratified weakly absorbing gas. We suppose a monochromatic light source is lying 
infinitely far away {z — — oo) and is producing a plane wave propagating in the +z direction. The needed properties 
of a stratified gas filling the spacetime are summarized by the density p(z) and the attenuation coefficient n u {z). The 
physical metric is 

ds 2 = -c 2 dt 2 + dx 2 + dy 2 + dz 2 , (42) 
and the vector fields k a and l a in Eq. ([7]) are trivially (u a = <5°/c) 

k a = k 0,0,1^ , 

r = ^(4<°'°4 (43) 

Equation (p~0|) then gives 

d~s 2 = -c 2 dt 2 + e 2a{z) (dx 2 + dy 2 ) + dz 2 , (44) 
where the conformal factor is related to the integral of the opacity by 

i r z 

a(z) = - a(z')dz'. (45) 
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We can now give a geometrical interpretation of Eq. (|35|) . Because the rays in Eq. (|43[1 are all running parallel to the 
z-axis the cross-sectional area of any bundle of rays is expanded by a factor e 2a ^ z ' for the conformally transformed 
spacetime, Eq. (|44|) . relative to initial the physical spacetime, Eq. (|42l) . From the definition of flux, we can conclude 
that it will be reduced by the reciprocal factor e _2<T . The reason the conformal factor was not applied to the whole 
4-D metric now becomes clear: a e 2<T factor in front of the entire physical metric would produce an undesirable time 
dilation factor e a as well as a frequency shift factor e _CT . These two combined would have contributed another factor 
of e~ 2cr to the wave's flux in the conformal spacetime. We can also look at the electric field E and the magnetic field 
H in both spacetimes and see that in the conformal spacetime they are each diminished by a factor e~ a . Consequently 
the flux, E x H, is reduced by e~ 2a . (This reduction in amplitude/intensity matches that of weak absorption given 
in but, not unexpectedly, differs from that of strong absorption, i.e.,Eqs. (28)-(30) of 0). 



V. ABSORPTION IN FLRW SPACETIMES 



The interesting example for cosmology is absorption by the Intergalactic Medium (IGM) which is naturally modeled 
as absorption by the cosmological fluid in Friedman-Lemaitre-Robertson- Walker (FLRW) spactimes. The familiar 
Robertson- Walker metric can be written in co- moving coordinates {u a = <5£/c) as 



ds 2 



-c 2 dt 2 + R 2 {t) 



dr 2 



1 — kr 2 



+ r 2 (d0 2 + sm 2 9d(/) 2 ) 



(46) 



where k = 1, 0,-1 for a closed, flat or open universe, respectively. The eikonal of Eq. (|32|) for a point source located 
at r = (see 12) is 



S(t,r)=R 



cdt' 
W) 



- sinn 1 [rl 



where 



!sin[r] k = +1 
r k = 

sinh[r] k = — 1, 

and Rq is the radius of the universe at the observation time to . The corresponding radial null geodesies are 



for which Eq. ([7]) gives 



(47) 



(48) 



(49) 



R(t) I 1 vT - kr 2 
~2Ro~ l~c' R(t) 



0,0 



and from which it follows that 

2t( a kb) = diag 
Equation (flU)) then gives the optical metric 



c 2 ; R 2 (t) 

' 1 — kr 2 



,0,0 



ds 2 



-c 2 dt 2 



R^t) | + e 2 °r 2 (d9 2 + sin 2 6 # 2 )| 



(50) 



(51) 



(52) 



An observer at (r, t) will see the source redshifted an amount z as a function of co-moving radius r and observing 
instant t, which can be obtained by integrating Eq. (|4"9"| . If the dynamics of a Robertson- Walker metric is determined 
by general relativity, and the gravity source is a mixture of non-interacting perfect fluids including a cosmological 
constant, the cosmology is called FLRW and has a redshift which can be found by inverting 



r(t, z) = sinn 





f z dz' ' 


lR(t)H(t) J 


o h(z')\ 



(53) 
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where 

h(z) = v/fi A + fi fe (l + z) 2 + fi m (l + z) 3 + fi r (l + z)±. (54) 

Here Q^,i} m ,n r are the respective density parameters for the assumed vacuum energy, cold matter and radiation 
content of the current universe. 

For spatially homogeneous and nondispersive (grey) absorption the differential optical depth dr is related to the 
opacity a by 

dr = —a(cdt) 

C 1 dz 

= a HjrTz)h(z-y (55) 

and the total optical depth r(r, t) of a source at (0, t e (r, t)) seen by an observer at (r, t) is the integral 

t{z) =w)L d+1w) rfz ' (56) 

The conformal factor a of Eq. ([52")) is therefore one-half this value. 

The presence of absorption changes the distance modulus-redshift relation n(z), e.g., the magnitudes of supernovae 
are corrected for absorption by the IGM before drawing a Hubble diagram. From the definition of distance modulus 

, = 5 log ^+ 25, (57) 

and luminosity distance g?l in terms of flux S 

/ L \l/2 

dL = [^s) < ^ 

we immediately obtain via Eq. (|39[) the absorption-corrected luminosity distance 

d L = e r/2 d L , (59) 

and the corrected distance modulus fl 

^^ + 2nho T - (60) 

This is the normal interpretation of the increase in luminosity distance caused by light absorption. If we use the 
optical metric to compute luminosity distance (see e.g. Q) we find that compared with the physical metric the 
cross-sectional area of the ray bundle is expanded by a factor e 2<T which reduces the received flux by a factor of e _2(T , 
i.e., from Eq. ([52"]) 

d L = e CT ( r -*°)(l + z)i? r = e a d L . (61) 

We see that by identifying a with r/2 we have the same luminosity distance in the optical spacetime with no absorption 
as in the physical spacetime with absorption. 

VI. INCORPORATING BOTH REFRACTION AND ABSORPTION INTO THE OPTICAL METRIC 

Up to now, we have been considering only absorption, however, we know that besides having its flux reduced, a 
light wave's path and speed can be altered (refracted) due to the presence of a polarizable material. An interesting 
and useful tool to study refraction in curved spacetime was developed by Gordon [l[ (see also Ehlers [§]). Gordon 
modified Einstein's physical spacetime metric to include the effects of a nondispersive refractive material on Maxwell's 
theory. The theory accounts for any polarizable material whose constitutive properties are summarized by two scalar 
functions, a permittivity e(x a ) and a permeability n(x a ). His optical metric g a b [l) is defined on the same differentiable 
manifold as Einstein's spacetime metric g a b and is related to it by 

gab = gab + (l - ^) u a u b , g ab = 9 ab + (1 - n 2 ) u a u\ (62) 
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where n(x a ) = ^fejl is the refractive index and u a is the 4- velocity of the optical fluid, normalized using the physical 
spacetime metric g a b- Null geodesies of the g a b metric are identical to light paths (timelike) in physical spacetime filled 
with a refracting material of index n{x). Even though we can incorporate absorption of only geometrical optics waves 
into Gordon's metric, his index of refraction theory applies to all Maxwell fields, see Q for some related examples. 
Gordon showed that solutions to Maxwell's theory (F a b, H ab ) in the presence of such a material can be found by 
solving a slightly modified set of vacuum (e = fi = 1) Maxwell equations for (F a b, F ab ) in his optical spacetime. The 
fields are connected by 

JJ a b — F ab 

A* 

F ab = F ab , (63) 

and the modified Maxwell equations are 2 



(64) 

We now generalize Gordon's optical metric to include effects of an additional isotropic and frequency independent 
(grey) attenuation coefficient k„ on a radiation field described by the geometrical optics approximation. To obtain 
this result we simply apply a transverse conformal transformation of Eq. (fT0|) to Gordon's metric Eq. (|62|) . The 
geometrical optics theory in Gordon's optical spacetime is almost identical to the geometrical optics theory of the 
physical spacetime given in Sec. lIIII (see Sec. Ill of @), e.g., the contravariant components of the Maxwell field are 
correctly given by Eq. (|32j) 




F ab = -2 Re{£ e lS/x °U a e b] }. (65) 

Here 

k a = g ab d b S, (66) 

is tangent to null geodesies of the optical metric (and tangent to the corresponding timelike "slower than light" 
curves of the physical metric), e b is the unit (in the optical metric) polarization vector, and S is the eikonal function. 
The significantly different equation is the transport equation for the amplitude of the wave given by Eq. (|3"Tj) , which 
becomes 

£ + 6£ + 4>S = 0, (67) 

(see Eq. (18) of [2|). The presence of the additional term <p = (1/4) ln(e//i) is due to the afore mentioned modification 
of Maxwell's equations [see Eq. (f6~4"|) ] in Gordon's optical spacetime. 
From Eqs. (fTO]) and (fTTj) the new optical metric becomes 

~9ab = e 2 °g a b + (1 - e 2a )2£ (a fc b) , g ab = ^g ah + (1 - e" 2 ^^. (68) 
where l a from Eq. (JT]) is a null vector field in Gordon's optical spacetime defined by 

k a u a 



2(u b k b ) 2 {u b k b Y 

and u a the fluid's 4-velocity normalized using Gordon's optical metric 



(69) 



u a - (70) 
n 



2 We thank John Ralston for pointing out to us that the modified equations can be derived from the Lagrangian C = 
-±y/-detgy/7/JlF ab F ab . 



10 



The reader should observe that the transverse conformal transformation does not alter the orthogonality of the two 
2-D subspaces nor does it alter the metric structure of the timelike 2-D space spanned by u a and k a . Rewriting the 
new optical metric Eq. (|68|) in terms of the physical metric g a b and physical observer u a we have 



Jab 



( 1 \ 1 - e 2<T 

e 2 °#ab + e 2<7 1 g ) U a U b + ~( Tx2 [ k ah + 2(u • k)k( a U b )] , (71) 



with inverse 



where 



and determinant 



The new Maxwell field 



~ab = £ -2a g ab + Q _ ^-2*^6 + 1 e * + 2n 2 (u ■ k)u {a ~k% (72) 

n z [u ■ ky 



k a = g ab S, b = 9 ab S,b = g ab S, b + (1 - n 2 )u a (u b S,b), (73) 



e 4<T 

det g = —5- det g. (74) 



pab _ ^ e -2a H ab^ 

F ab = F ab , (75) 

satisfies the same equations as F ab , i.e., Eq. ([64]) . but has the advantage of incorporating both refraction and 
absorption. We now have a correspondence of geometrical optics waves in two spacetimes, the physical and the 
optical. In the physical spacetime the wave travels at a reduced speed c/n with an intensity that is reduced by 
absorption as in Eq. (|37| whereas in the optical spacetime the wave travels at speed c with no extinction. 



VII. THE OPTICAL METRIC IN FLRW COSMOLOGY WITH BOTH REFRACTION AND 

ABSORPTION 

The formalism developed above appears complicated when applied to an arbitrary spacetime, however, for specific 
cases the formalism is more transparent. For example we reconsider FLRW spacetimes, but this time with both 
refraction n and absorption k associated with the cosmic fluid. The result is elegantly simple. Gordon's optical metric 
for the Robertson- Walker metric 0] is 



c 



2 



dr 



ds 2 - = -—dt 2 + R 2 (t) \ 1 _ kr2 + r 2 {d6 2 + sin 2 6 dc/) 2 ) } . (76) 



This metric does not measure distance and time, that is reserved for Eq. (J46J). A dynamical theory that produced a 
gravitational field containing the index of refraction as given in Eq. (|76p would be quite strange in that it would be 
partially sourced by electromagnetic polarization densities; which is not stress, energy, or momentum. The significant 
property of the optical metric that we use here is that its null geodesies are identical with the speed c/n photons of 
Eq.gBJ. 

For spherical waves the radial null geodesies are 

k a = Ro[ ,0,0 1, (77) 

and from Eq. (|69|) we find 

i a = Hr\ --, v ^ ,o,ol ■ ("si 





Immediately we find 
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which gives us the new optical metric 



dS 2 = -^dt 2 + R 2 (t) |y^2 + e 2a r 2 {d9 2 + sin 2 6 dtf) ) . ( SO) 



n 

With dynamics supplied by the FLRW solutions the scalar function a in the above is 

cdf 



n 

1 C a(z') , 

~2H(f)J n(z>){\ + zi)h{z>) ( 1} 

Note the difference between wavelength redshift 2 and frequency rcdshift z„ (see 0): 

1+z = wry 1+Zn = ^tT)W (82) 

We have discussed the impact of light refraction on the distance- redshift relation in [2] ■ Now incorporating absorp- 
tion we find from Eq. (|80p that the apparent size distance d A to the source in the new optical metric is 



d A = R{t e )e<°' u h = R(t e )r = d A , (83) 
is the same as the apparent size distance in Gordon's spacctimc 



7 / x 1 C 1 . 

d A (z n ) = — — - — — — 7= Sinn 
[1 + z(z n )\ H y/p^ 



z{z ^ dz 



n(z)h(z) 



(84) 



Finally the luminosity distance d~L is (see 

d L = (1 + z n ){\ + z)e a d A = e T l 2 d L . (85) 

We have thus obtained the same luminosity-redshift relation we obtained in for the real eikonal case (refer to 
Eq. (66) of 0) by changing the geometry rather than absorbing some of the wave's intensity. Just as in 0, Eq. l[55|) 
differs from the standard reciprocity relation 0, [l(J fTl| 

d L {z) = (l + z) 2 d A (z), (86) 

in two aspects: first, the existence of refraction causes photon orbits to deviate from null geodesies in the physical 
spacetime, this giving the 1 + z n factor instead of 1 -I- z; second, light absorption (expressed by the conformal factor 
a) violates the photon number conservation law. For an interesting discussion about using reciprocity relation as a 
probe of acceleration/exotic physics, see [l2j . 



VIII. AN APPLICATION: FLRW SPACETIME WITH BOTH REFRACTION AND ABSORPTION 

In 0, [H[ we have interpreted the observed apparent increase in the universe's expansion rate as caused by light 
refraction and absorption respectively, instead of by a cosmological constant A. In [131 ] we used a two parameterpure 
refraction model, i.e.,n(z) = 1 + pz 2 + qz 3 , without absorption to fit the supernovae gold sample In we 

fit the sample with a pure absorption model, i.e., we took n = 1, and the opacity a(z) = p(z)n(z) = const. In this 
section, we fit the same data set [3, [TH, [IH, [T3, 13 with a cold dark matter model containing both refraction and 
absorption. We use simple expressions n{z) = 1 + pz 2 and a(z) — const that depend on only two parameters. Since 
we are concerned with the matter dominated era, we have excluded radiation (Cl r — 0) and since we are trying to 
only emulate acceleration we take A = 0. We choose the current Hubble constant to be Hq = 65km/s/Mpc and fit 
our two parameter (a, p) model for different choices of f2 m . The scaled Hubble function h{z) in this case simplifies to 

h{z) = (l + z)y/l + Q m z. (87) 

The refraction and absorption corrected distance-redshift relation is now written as 

c e <K*») . H 2 ") dz' 

d L {z n ) = (1 + z n ) — -= S mn^/\Sh\ / - (88) 

H Jo n[z')h{z') 
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where 



and 



a c f z{Zn) dz' 
a{Zn) =2H~ J {l + z')n{z>)h{z>y m 



l + z n = -±±^. (90) 
1 + pz* 

We use the 178 supernovae from the gold sample [l9[ with rcdshifts greater than cz = 7000 km/s to avoid any 
Hubble bubble. Our results are shown in Figures [T] and In Fig.[TJ we show the Aju(z) curves for different model 
parameters. Here A/i(z) = /i(z) — p,p(z), where /j>f{z) is the distance modulus of the fiducial, dark matter only model 
(horizontal black dashed curve in Fig.[T]), i.e.,£lA = 0, fl m — 0.3, a = 0,p = 0. In each of the four frames the green 
curve (grey in black & white) is the concordance model, fl\ = 0.7, fi m = 0.3, a = 0,p = 0. In the upper left panel 
we show absorption models with no refraction index (Q m < 1 for all curves in this panel). The dotted red curve 
has: £1a = 0, Q m = 0.05, a = 0.7 x 10~ 4 Mpc -1 , n = 1. The short-dashed red curve has: = 0, Q, m = 0.3, a = 
1.3 x 10- 4 Mpc~\n = 1. The solid black curve has: tt A = 0,fi m = 0.73, a = 2.2 x 10~ 4 Mpc _1 , n = 1. In the upper 
right panel the models are flat. The blue curve (bottom curve) has: Q\ = 0, Q m = 1.0, a = 2.5 x 10 _4 Mpc ,p = 0. 
The red curve (top curve) has: — 0, Q m — 1.0, a = 2.5 x 10 _4 Mpc~ ,p = 0.024. In the bottom left panel the 
blue curve (bottom curve) has: fl\ — 0, O m = 1.5, a = 3.3 x 10~ 4 Mpc -1 ,p = 0. The red curve (top curve) has: 
f?A = 0, r2 m = 1.5, a = 3.3 x 10~ 4 Mpc -1 ,p = 0.042. In the bottom right panel the blue curve (bottom curve) 
has: Q\ = 0,O m = 2.0, a = 4.1 x 10 _4 Mpc ,p = 0. The red (top curve) curve has: Q\ = 0, Q m = 2.0, a = 
4.1 x 10- 4 Mpc" 1 ,p = 0.049. 

In Fig. [2] we show the confidence contours (68.3%, 95.4%, and 99.73%) of our two parameter (a,p) model with 
different choices of fi m . We show 6 different cases with f2 m = 0.05 (baryonic matter only), f2 m = 0.3 (dark matter 
only), £l m = 0.73 (our best least \ 2 P ure absorption model, see Q), fi m = 1.0 (flat universe), Q m = 1.5 and 2.0 
(closed models). We restrict the parameter p to be nonnegative to keep the light speed v — cjn less than c. For 
the first three cases, i.e., Sl m = 0.05,0.3 and 0.73, the best fits occur for p = Q, which suggests that introducing 
additional refraction parameters would not significantly improve the fitting when absorption is present. However, as 
f2 m increases, nonvanishing p values do give better fits. For f2 m = 1.0, the best fitting parameters (x 2 = 1.04) are 
a = 2.5 x 10" 4 Mpc _1 ,p = 0.024. For fi m = 1.5, the best fitting (x 2 = 1.05) parameters are a = 3.3 x 10" 4 Mpc _1 ,p = 
0.042. For Q m = 2.0, the best fitting (x 2 = 1.06) parameters are a = 4.1 X 10~ 4 Mpc -1 ,p = 0.049. This can also 
be seen from Fig.[TJ In the upper right and two bottom frames, the blue and red curves have the same absorption 
coefficient a, the difference is that the red curve has p nonzero whereas the blue curve has p = 0. The inclusion of p 
for these large f2 m cases improves the fits. 

The Supernova data is currently considered to be the most compelling evidence for the existence of dark energy 
because of the apparent acceleration observed in the expansion of the universe (see e.g. I20I [pill') Competing 
interpretations have been proposed, e.g ., evolutionary effects [H, [23|, local Hubble bubbles [24 l25l | , a bsorption 
[H, [U [H[ , modified gravity [2^, [3(| [Hf and others such as slowly changing fundamental constants [U Hi| . 

For supernovae there are at least four different sources of opacity; the Milky Way, the hosting galaxy, intervening 
galaxies, and the IGM that should be taken into account. The Galactic absorption has been studied extensively 
[34l [HI and early constraints on properties of IGM were often obtained assuming the applicability of Galactic dust 
properties (see e.g. [HJ). The luminosity of high redshift supernovae have been corrected for host galaxy absorption 
using the Galactic reddening law, see e.g., [13, HH, H3, HI] . As has been pointed out caution should be exercised when 
applying Galactic dust properties to the IGM, since the composition, size, shape, and alignment of intergalactic dust 
could be significantly different than that of Milky Way dust. Aguirre [2y, [27| introduced a carbon needle model 
and showed that dust grains of larger size (> 1 /xm), which should be preferentially ejected by star burst galaxies, 
would have relatively higher opacities (e.g. k ~ 10 5 cm 2 j _1 ) and much greyer absorption curves, and therefore 
mig ht escape the reddening censorship based on Galactic reddening law. Evidence against this model appeared in 
[39l |40|. l4l| . Goobar et al. [12] introduced a replenishing dust model in which the dilution caused by cosmic expansion 
is continually replenished (p = po(l + z) 3 for z < 0.5, p = po for z > 0.5). This dust model is indistinguishable from 
Oa and cannot be ruled out by Supernovae data alone (see table 5 of [lB]). Bassett & Kunz [l§] claimed to rule out 
the replenishing model at more than 4<r by considering violations of distance duality (reciprocity relation) . Ostman & 
Mortsell [43| further constrained the magnitude of grey dust absorption from Quasar colors and spectra and claimed 
that for a wide range of intergalactic dust models, extinction larger than 0.2 mag is ruled out (see also 0, EH). A 
new grey dust model is proposed in [46j. More complicated and fine tuned dust models will keep emerging in the 
future until dark matter/energy has been identified, if in fact it exists. 



13 




0.02 0.050.100.20 0.501.00 




0.02 0.050.100.20 0.501.00 




0.02 0.050.100.20 0.501.00 




0.02 0.050.100.20 0.501.00 



FIG. 1: A/x versus z. In each of the four frames the fiducial model (horizontal black dashed) is the now disfavored dark matter 
only model, i.e., Qa = 0, O m = 0.3, a — 0,p = 0, and the green curve (gray in black & white) is the concordance model, £7a = 
0.7, fi m = 0.3, a — 0,p = 0. The upper left panel contains pure absorption models with no refraction: the dotted red curve has 
fi A = 0, tom = 0.05, a = 0.7 x 10~ 4 Mpc _1 ,n = 1; the short-dashed red curve has fi A = 0, fi m = 0.3, a = 1.3 x 10" 4 Mpc"\ n = 1; 
the solid black curve has Qa = 0, fi m = 0.73, a = 2.2 x 10 _4 Mpc _1 , n = 1. In the upper right panel: the blue (bottom) curve 
has Q A = 0, £l m = 1.0, a = 2.5 x 10~ 4 Mpc"\p = 0; the red (top) curve has fi A = 0, fi m = 1.0, a = 2.5 x 10~ 4 Mpc _1 ,p = 0.024. 
In the bottom left panel: blue (bottom) curve has fi A = 0, Cl m = 1.5, a = 3.3 x 10~ 4 Mpc~ 1 ,p = 0; The red (top) curve has 
Qa = 0, fi m = 1.0, a — 3.3 x 10 _4 Mpc _1 ,p = 0.042. In the bottom right panel: blue (bottom) curve has Qa = 0, fi m = 2.0, a = 
4.1 x 10 _4 Mpc _1 ,p = 0; the red curve has Q A = 0, Q m = 2.0, a = 4.1 X 10~ 4 Mpc _1 ,p = 0.049. 



IX. DISCUSSION 



In this paper, we introduced the concept of a "transverse" conformal transformation, which allowed us to equate 
the intensity reduction of light caused by absorption with a geometrical reduction caused by the expansion of the 
cross-sectional area of light ray bundles without changing other desirable properties of waves, such as frequency and 
wavelength. This application of conformal transformations is new. We used it to generalize Gordon's optical metric 
to include light absorption via such a conformal transformation. This generalization is fundamentally different from 
in which we included light absorption by making Gordon's metric complex. In the complex metric formalism, we 
distinguished two different cases: strong and weak absorption. The strong case is where a non-negligible amount of 
absorption occurs on a wavelength scale and the weak case is where absorption is significant only over multitudes of 
wavelengths. In the strong case the eikonal remains complex but in the weak case it can be taken as real. In this 
paper we were able to replace the effects of weak absorption by a special conformal transformation. The disadvantage 
of this formalism is that the optical metric depends on the eikonal of the wave itself whereas in [2| it did not. We 
used this new optical metric to derive the absorption and refraction corrected distance-redshift relation for FLRW 
spacetimes and obtained the same expression as in Q for weak absorption. In [l3| we fit supernovae data with a pure 
refraction model, and in we fit it with a pure absorption model. In Sec. lVIHl we fit this data with a cosmological 
model possessing both refraction and absorption. We have shown that a single parameter polynomial approximation 
for the refraction index, i.e.,n(z) = 1 +pz 2 , together with a constant opacity parameter a fits the data well. More 
realistic z dependent models for a and n would be appropriate. We have not proposed a physical source for the 
needed refraction and/or absorption. What we have done here is to construct a geometry modification whose effects 
are equivalent to absorption. It is not a new absorption theory only a new way of looking at its effects. It equates the 
decrease in intensity of a wave caused by absorption with a decrease caused by a change in the underlying spacetimc, 
i.e., a transverse conformal transformation. We used this modified geometry to compute the distance-redshift relations 
in FLRW and applied it to the Hubble curve for type SNe-Ia. In this application we ignored effects of inhomogeneities 
including the possibility of having a multi-valued conformal factor. We assumed that none of the SN la used were 
strongly lensed. Since the flat concordance model is consistent with other observations, e.g., the angular position of 
the first acoustic peak as measured by WMAP (Wilkinson Microwave Anisotropy Probe) team [4^, [48| , and baryonic 
acoustic oscillations (BAO) detected in galaxy surveys [4^, [5(| [EH , an absorption and/or refraction theory cannot be 
on firm ground unless it is also consistent with these additional observations. We leave these and other applications 
to future efforts. 
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FIG. 2: Confidence contours (68.3%, 95.4%, and 99.73%) for each fixed f2 m model. The x coordinate is the absorption coefficient 
a in unit of 10 _4 Mpc _1 , and the j/-axis is the refraction index parameter p (n = 1 + pz 2 ). In the top left panel, fi m = 0.05, 
with least chi-square (per degree of freedom) x 2 = 1-09, the best fitting parameters are a = 0.7 x 10 _4 Mpc _1 , p = 0. In the 
top right panel, f2 m = 0.3, least x 2 = 1-05, the best fitting parameters are a — 1.3 x 10 _4 Mpc _1 , p = 0. In the middle left 
panel, £7 m = 0.73, least x 2 = 1-035, the best fitting parameters are a = 2.2 x 10 _4 Mpc _1 , p — 0. In the middle right panel, 
n m = 1.0, least x 2 = 1-042, at a = 2.5 x 10~ 4 Mpc~\ p = 0.024. In the bottom left panel, fi m = 1.5, least x 2 = 1-05, the best 
fitting parameters are a — 3.3 x 10 _4 Mpc _1 , p = 0.042. I the bottom right panel, Q m = 2.0, least x 2 = 1-06, the best fitting 
parameters are a = 4.1 x 10 _4 Mpc _1 , p — 0.049. 
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